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STIMULATED  RAMAN  SCATTERING  OF  INTENSE  LASER 
RADIATION  IN  A  NUCLEAR-DISTURBED  SPACE  PLASMA 

I.  INTRODUCTION 

A  key  component  of  proposed  Strategic  Defense  Initiative  (SDI)  systems 
is  high  pover  coherent  radiation  sources,  either  ground  or  space-based. 
The  propagation  characteristics  of  such  radiation  beams  in  the  ambient  and 
nuclear-disturbed  environment  is  of  the  utmost  importance  for  the 
successful  implementation  of  SDI  goals  and  objectives. 

Much  work  has  been  devoted  to  the  study  of  laser  propagation  in  the 
ambient  atmosphere^ *^.  For  low  beam  pover  densities,  veil  knovn  effects 
such  as  turbulence  broadening  of  the  beam,  distortion  of  beam  spatial 
coherence,  and  random  refraction  and  vander  of  the  beam  have  been  studied 
both  experimentally  and  theoretically.  At  higher  beam  pover  densities, 
nonlinear  effects,  e.g.,  thermal  blooming,  stimulated  scattering,  and 
optical  breakdovn  can  play  an  important  role  in  beam  propagation  by 
enhancing  beam  absorption  and  scattering. 

For  the  nuclear-disturbed  environment,  vhile  many  studies  have  been 
made  to  assess  the  impact  of  the  high  altitude  nuclear  explosion  (HANE) 
environment  on  the  propagation  of  radiation  in  the  radio  frequency  range, ’ 
very  little  vork  has  been  performed,  particularly  at  high  beam  pover 
densities,  dealing  vith  laser  and  optical  propagation  through  the  HANE 
environment.^-^  In  this  study  ve  focus  on  plasma  parametric  instabilities 
driven  by  the  radiation  beam.  These  instabilities  are  dependent  on  the 
beam  pover  density  and  can  affect  beam  propagation  by  scattering  beam 
radiation  energy.  The  fastest  groving  plasma  parametric  instabilities  will 
involve  high  frequency  electron  dynamics.^  From  the  results  of  our  initial 
vork^  ve  focus  on  stimulated  Raman  scattering  (SRS)  since  this  instability 
can  be  quite  rapid  and  can  grov  on  fast  «pe-1  tine  scales  vhere  u>pe  is  the 
electron  plasma  frequency.  The  outline  of  this  report  is  as  follovs.  In 

Manuscript  approved  September  22,  1988. 
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Section  II  ve  outline  the  theory  and  attendant  assumptions  in  stimulated 
Raman  scattering  in  nuclear-disturbed  plasmas.  In  Section  III  ve  present 
our  results  giving  growth  rates  of  stimulated  Raman  scattering  (SRS)  for  a 
typical  single  HANE.  In  addition  ve  also  compute  growth  rates  for  SRS  for 
the  multi-burst  case.  Finally,  in  Section  IV  ve  summarize  and  discuss  our 
results. 

II.  THEORETICAL  MODEL 

The  basic  geometry  of  our  model  consists  of  a  monochromatic 
electromagnetic  vave  vith  frequency  «0  and  wavenumber  kQ  propagating  in  a 
plasma  imbedded  in  a  magnetic  field  as  shown  in  Fig.  1.  Since  ve  are 
considering  frequencies  Uq  >>  «pe,  2e  the  effects  of  the  ambient  magnetic 
field  on  the  electromagnetic  wave  will  be  small;  Here  Wpe  is  the  electron 
plasma  frequency  and  2g  the  electron  gyro frequency.  In  addition,  since  ve 
will  be  investigating  the  scattering  of  high  frequency  modes  vith 
frequencies  |o>|  1  «pe  »  2e  »  »  2^,  i.e.,  Langmuir  waves,  ve  take,  as 
a  first  approximation,  the  electrons  and  ions  to  be  unmagnetized.  Here  copj 
is  the  ion  plasma  frequency  and  2^  the  ion  gyrofrequency.  As  a  result,  the 
general  equilibrium  implied  by  our  model  consists  of  electrons  oscillating 
with  high  velocity  vQ  -  eE0/rae»0  in  a  stationary  ion  background  where  e  is 
the  electron  charge,  E0  the  electric  field  associated  vith  the  incident 
electromagnetic  vave,  and  me  is  the  electron  mass. 

If  one  perturbs  this  equilibrium  vith  a  density  and  electric  field 
perturbation  of  the  form  exp  [i(k^  •  x  -  w^t)]  then  beat  oscillations  will 
form  with  frequency  «2  =  ±  %  and  vave  number  |S2  =  iSl  *  iSo *  These 
sidebands  will  in  turn  beat  vith  (u>0,  kQ)  and  cause  the  growth  or  decay  of 
the  mode  (ttj,  k^)  provided  the  Manley-Rove  relations 
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«o  =  +  o»2  (la) 

£0  *  +  l£2  Ob) 

are  satisfied.  As  a  result,  a  Langmuir  vave  and  a  sideband  electromagnetic 
vave  can  be  parametrically  driven  to  large  amplitudes  at  the  expense  of  the 
pump  electromagnetic  vave. 

When  (!>}■<»  and  k^  ■  k  are  Langmuir  waves,  the  dispersion  relation  for 
stimulated  Raman  scattering  can  be  found  using  the  combined  Vlasov-Maxvell 
equations  and  standard  techniques  to  yield^*® 

X~  +  xprr  *(2kovosin2<x  c°s2e)/»0(w  -  Att  +  iYd)  (2) 

where  a  is  the  angle  between  the  beat  wave  k  -  k0  and  vQ  (see  Fig.  2),  9  is 
the  angle  between  k  and  kQ,  Cm  ■  (c2/«0)(k’k0  -  k2/2),  y<j  is  the  damping 

rate  of  the  free  electromagnetic  wave  with  y<j  *  ve“pe/2wo'  ve  is  the 
electron  collision  frequency,  and  Xe(Xi>  is  the  electron  (ion) 
susceptibility  given  by 

Xe(o»,k)  -  2  (k^)"2 [l  +  CeZ(Ce)][l  +  (ive/kve)Z(Ce)]~1 

and 
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X^k)  -  2(Te/T.)<kXD)'2[l  +  qZCCj)] 


with  £g  -  («  +  ive)/kve,  Ci  ■  w/kv^,  Xg  is  the  Debye  length,  ve(vj)  is  the 
electron  thermal  velocity,  Z  is  the  plasma  dispersion  function  and  Te(T^) 
is  the  electron  (ion)  temperature.  Ve  note  that  the  Manley-Rove  relation 
requires  k  »  2kocos0  in  Eq.  (2). 

Bq.  (2)  has  been  solved  for  u  >  ^  ^  iy  in  several  limits.  For 
example,  if  kXg  »  2kocos0Xg  «  1,  so  that  the  plasma  wave  is  not  heavily 
damped,  it  is  found2'®  that  Mj.  *  (“pe  ♦  k2v|)^2,  the  Bohm-Gr oss  frequency, 
and 


r  ■ 


-  hrP 


V 


I  [<■ 


4(vo/c)2sin2<x  cos20  w0wpe]1/2  (3) 


where  Yp  is  the  total  damping  (Landau  and  collisional)  rate  of  the  free 
plasma  wave.  The  threshold  beam  power  density  is 


-2  -2 

sin  a  cos  QCYp/WpgHrj/^) 


(4) 


giving  maximum  growth  rate  for  rmax  far  above  threshold  (yp  <<  y  <<  Mpe) 


Ymax  =  (v0/c)|sinajcos9(w0a)pe)1/2 


(5) 


In  this  case  the  grovth  rate  assumes  a  linear  dependence  on  E0  or  v0  if 
(vQ/c)  >  (ypIsinalcosSfcipgttQ)!^  j .  On  the  other  hand,  for  short 
wavelengths  kXg  -  2kocos0  »  1,  is  can  be  shown?’®  that 

«  »  +  2(vQ/c)2sin2a  cos20  “0 [xg(&«»2kocos0)/(l  +  Xe(A»,2kocos9)) j 

giving  a  grovth  rate 
Y  *  2(vQ/c)2sin2a  cos20  »oIm |xe( Aw,  2kocos0)/(l  +  Xe(&»t  2kQcos0))J-  rd 

(7) 

If  collisional  damping  of  the  electrostatic  wave  is  neglected,  Eq.  (6)  and 
(7)  give  Wj.  -  kve  and 


Y  -  («0/2)(vo/c)2(n/2)1/2sin2«  exp(-  l/2)(kQXd)-2  -  Yd  (8) 

with  threshold 

<Vo/c)TH  *  (8/H)1/2(Yd/«0)exp(l/2)(koXD)2sin'2a  (9) 

and  maximum  grovth  rate 

Ymax  ’  («o/2)(v0/c)2(K/2)1/2sin2a  expC-l^Mk^)-2  (10) 

For  the  nuclear  environment  the  quantity  can  be  of  order  unity  with 

the  result  that  Eq.  (2)  must  be  solved  numerically. 
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III.  RESULTS 


We  now  proceed  to  solve  Eq.  (2)  using  parameters  typical  of  HANE 
plasmas.  For  a  CV  laser  with  108V  peak  power  with  a  5-10  m  beam  diameter, 
we  find  beam  power  densities  I  =  102  V/cm2  and  resulting  quiver  velocities 
v0  •  25.6  [I(V/cm2) ]*72  X^Cpm)  cm/sec  =  2.5  x  103  cm/sec  where  X^  is  the 
laser  wavelength.  This  gives  nonrelativistic  quiver  velocities  vQ/c  a  10“7. 
For  the  pulsed  mode  with  108J  of  energy  with  repetition  rate  of  5u  sec,  we 
find  v0/c  a  10"5. 

Fig  3(a)-(b)  gives  the  growth  rate  y  for  stimulated  Raman  scattering  as 
a  function  of  vQ/c  for  both  the  pulsed  and  CV  modes  using  parameters 
typical  of  the  later  time  (t  30  min)  evolution  of  a  single  high  altitude 
nuclear  burst.  At  late  times,  the  HANE  plasma  can  extend^  several 
thousands  (hundreds)  of  kilometers  parallel  (perpendicular)  to  the 
geomagnetic  field  and,  as  a  result,  provide  a  large  plasma  volume  for 
scattering  to  occur.  We  take  ne  *  107  cm-3,  Te  =*  Ti  *  0.2  eV,  and  consider 
both  1  ym  (kQXj)  -  6.2  x  103)  and  10  urn  (^Xg  *  6.2  x  102)  laser 
wavelengths,  we  find  the  shortest  growth  times  y~*  =  1-10  sec  for  the  10 
pm  case  in  the  pulsed  mode. 

Fig.  4(a)-(b)  displays  the  growth  rate  for  stimulated  Raman  scattering 
as  a  function  of  vQ/c  again  for  both  the  pulsed  and  CV  cases  but  using 
parameters  typical  of  computer  simulations  [R.  Dana,  private  communication! 
of  a  multiple  high  altitude  nuclear  burst  scenario.  Here  we  take  ne  =  2  x 
1011  cm-3,  Te  »  'c  eV,  and  consider  both  the  10  urn  (kQXg  =  10)  and  1  um 
(k0Xg  »  100)  cases.  Here  with  find  much  larger  growth  rates  with  the 
fastest  occurring  on  time  scales  y_1  =  10“  3  -10-2  sec  for  the  pulsed  mode. 
The  growth  rates  y  and  corresponding  growth  times  for  the  CV  case  are 
smaller  due  to  the  reduction  in  beam  power  density. 
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IV.  SUMMARY  AND  DISCUSSION 

We  have  computed  the  growth  rates  y  of  stimulated  Raman  scattering  of 
both  10  um  and  1  urn  laser  radiation  using  plasma  parameters  typical  of  a 
single  and  multiple  high  altitude  nuclear  explosion.  Ve  have  considered 
both  CV  and  pulsed  modes  of  laser  operation  using  parameters  representative 
of  proposed  SDI  systems. 

Por  the  single  burst  case,  ve  find  relatively  long  e-folding  growth 
times  *  1-10  sec  in  the  late  time  regimes  for  both  the  10  um  and  1  pm 
cases  in  the  pulsed  mode.  The  decreased  power  densities  in  the  CV  case 
result  in  much  smaller  grovth  rates  and  longer  growth  times. 

On  the  other  hand,  for  the  multiple  burst  scenario,  where  plasma 
densities  are  expected  to  be  much  higher  than  the  single  burst  case,  we 
find  faster  growth  times  y~^  -  10~^  -  10“^  sec.  For  the  pulsed  10  um  case 
typical  e-folding  growth  lengths  Le  =  cy-1  are  on  the  order  of  100  km.  For 
nuclear  plume  dimensions  on  the  order  of  several  thousands  of  kilometers, 
one  would  expect  many  e-folding  lengths  over  which  stimulated  Raman 
scattering  can  grow.  Depending  on  the  angle  of  arrival  of  the  beam  with 
respect  to  the  HANE  plasma,  which  is  expected  to  extend  several  thousands 
of  kilometers  in  altitude,  Raman  scattering  of  beam  energy  in  the 
multiburst  cas*  is  expected  to  be  strong  with  subsequent  deleterious  beam 
attenuation,  erosion,  and  degraded  propagation. 

In  order  to  more  fully  assess  the  impact  of  plasma  parametric 
instabilities  on  beam  propagation,  several  remaining  issues  need  to  be 
addressed.  Ve  have  made  the  assumption  of  a  homogeneous  plasma.  In  all 
probability,  the  HANE  plasma  will  be  highly  structured  or  turbulent  on 
scale  sizes  comparable  to  the  beam  radius  and  smaller.  This  density 
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structure  will  affect  the  growth  and  evolution  of  Langmuir  waves  driven 


unstable  by  stimulated  Raman  scattering.  In  addition,  we  have  not 
discussed  nonlinear  saturation  mechanisms  for  stimulated  Raman  scattering 
in  a  nuclear  environment.  Resolution  of  these  problems  will  give 
definitive  information  on  beam  propagation  in  the  SANE  environment. 
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Fig.  1  Basic  geometry  of  electromagnetic  wave  with  wavenumber  k0,  wave 
electric  field  E0,  and  wave  magnetic  field  B0>  propagating  in 
plasma  embedded  in  magnetic  field  BgxT' 
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Fig.  2  Definition  of  angles  a  and  0  used  in  Eq.  (2). 
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RAMAN  SCATTERING 

Ne=lxl07cm"3 

PULSED 


A:Xl=  lA*m 

B:Al=  lCW* 


V„/C 

Plot  of  growth  rates  r/<*)pe  vs.  vQ/c  for  stimulated  Raman 
scattering  for  Ne  -  107  cm-3,  Te  *  0.2  eV,  ve  »  5  x  103  sec"1  with 
(a)  10"5  <  v0/c  <  10"4  and  (b)  10"7  <  vQ/c  <  10"6. 
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RAMAN  SCATTERING 
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Plot  of  growth  rates  r/«pe  vs.  v0/c  for  Raman  scattering  for  Ne 


2  x  1011  cnT^,  Te  =  2  eV,  \>e  =» 

IQ*4  and  (b)  10"  7  <  v./c  <  10"6. 
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Fig.  4.  (Continued)  Plot  of  growth  rates  y/«pe  vs.  vQ/c  for  Raman  scattering  for  Ne  * 

2  x  10*1  cm-3,  Te  =»  2  eV,  \>e  *  10^  sec-1  with  (a)  10~3  <  vQ/c  < 
10-4  and  (b)  1CT7  <  vn/c  <  10"6. 


